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Heat Removal from 
Nuclear Reactors 


?! 


For a reactor to operate in a steady state, with an internal temperature distribution 
that is independent of time, all of the heat released in the system must be removed 
as fast as it is produced. This is accomplished, in all reactors except those operating 
at very low power levels, by passing a liquid or gaseous coolant through the core 
and other regions where heat is generated. The nature and operation of this coolant 
system is one of the most important considerations in the design of a nuclear reac- 
tor. 

The temperature in an operating reactor varies from point to point within the 
system. As a consequence, there is always one fuel rod, usually one of the rods 
near the center of the reactor, that at some point along its length is hotter than 
all the rest. TMiffiMimum-fuete 

theireactf^^ of^Ae^uefeHowever, 
metallurgical considerations place an upper limit on the temperature to which a fuel 
rod can safely be raised. Above this temperature, there is a danger that the fuel may 
melt, which can lead to the rupture of the cladding and release of fission products. 
Qngffiij^^ a reactor coolant system is;to provide/ 

fo||^ desired power level while '^suring that 

tBe|maxim^ below this predetermined yalue. 
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It should be noted that, from a strictly nuclear standpoint, there is no theoretic 
cal upper limit to the power level that can be attained by any critical reactor having 
sufficient excess reactivity to overcome its negative temperature coefficient. This; 
by removing control rods and plaang a reactor on a positive period, its power could 
be increased indefinitely were it not for the fact that eventually a point would he 
reached where the coolant system is no longer able to remove all of the heat being ^ 
produced. Beyond this point, the fuel would heat up and eventually a portion of the^jj* 
core would melt down. This situation is avoided in the actual operation of a react^H 
by reinserting some of the control rods or adding boron to the shim system to return 
the reactor to critical when a desired power level has been reached. Throughout the 
present chapter, it is assumed that the reactors under cbhsidemdqn are critical and 
operating at a constant powers | | 

Before beginning the discussion of the design of reactor coolant systeiSE 
word is in order concerning units. As pointed out in Appendix I, where the matter/' 
of units is considered in some detail, the SI system of units has been adopted by, | 
most of the nations of the world, with the exception of the United States. All U.S.fff- 
scientists and many US. engineers and engineering societies also use this system. 
However, there is still strong resistance to the adoption of SI units. As a resul' 
most American manufacturers of nuclear power equipment base their designs on 
the English system. Since this state of affairs is likely to persist for several ye r 
to come, English units are used throughout the present chapter as the primary unitg||li 
SI units are provided where appropriate. Tables for conversion from English to 
SI units are given in Appendix I. 


8.1 GENERAL THERMODYNAMIC CONSIDERATIONS 


Ill 


From a thermodynamic point of view, a nuclear reactor is a device in which energy 
is'pntMucedrandL as indicated in Fig. 8.1, heat is 

released in a reactor at the rate of q BTU/hr or watts and absorbed by the coolant,; 
which enters the reactor at the temperature r in and exits from the reactor at the 
temperature T mu passing through the system at the rate "of w Ib/hr or kg/hr. | 

With all power reactors except the BWR and RBMK, there is no (net) chang 
in phase of the coolant as it passes through the reactor; that is, the coolant does nc 
boil. In these reactors, the heat from the reactor merely increases the temperature oh 
the coolant, a process that occurs at essentially constant pressure — namely, the re- 
actor coolant pressure. The "Heat in BTUs ^ or joules required tb;raise^ 


out. 
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Figure 8.1 Schematic drawing of 
coolant flow through a reactor. 


where c p (T) is the specific heat at constant pressure per unit mass of the coolant. 
Since w lb or kg of coolant flows through the reactor per hour, th^t^lwffitff 
heat;pfpduced; in the reactor is absorbed by r thexoolant is given by 


(8.1) 


Equation (8.1) can also be written in terms of the thermodynamic function 
enthalpy. This is defined for a substance by the relation 


(8.2) 


wfea^&sM^^Py per-unit mass of the substance (pi^tfcmthalpyfru<isf> 
fe&S&PlbS&^^'H^ mass,. P ,is the pressure,, and^is .the* specific, volume • 
of the substance (i,e:f me volume -per unit mass ft 3 /lb or m 3 /Eg). The units of h 
are Btu/lb or joules/kg. In thermodynamics, it is shown that whJlheafSadded to 
a sub^tai^eiltleonstant. pressure? essentially all-ofithe-heat is us6d to increase its 
fM^PK^ Thus, if and h oal are the specific enthalpies of the coolant entering 
and leaving the reactor, respectively, it follows that , / 

^ b(A+ foi** (8.3) 


1 A very small amount of heat may also be used to increase the kinetic energy of the substance 
or change its gravitational potential energy by virtue of raising the substance from a lower to a higher 
level in the earth's gravitational field. However, both of these contributions are negligible for reactor 
coolants. 
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It also follows that Eq. (8. 1 ) can be written as 


chap, ay 


$r w S h ™ r Afi). (8.4) 

The situation is somewhat more complicated in the case of a BWR or an 
RBMK. In these reactors, a portion of the water passing through the core is va-if! 
porized to steam, exits from the reactor via a steam pipe to the turbine, and laterfij 
returns as feed water from the condenser and reheaters. Most of the water that 
goes through the core is recirculated within the reactor. However, in steady-state 
operation, there is no net absorption of heat by the recirculating water. Therefore, 
although the incoming feed water mixes with the recirculating water before pass! 
ing through the core, the overall effect of the reactor is simply to vaporize the feedl 
water, as is indicated schematically in Fig. 8.2. f% 

The change in enthalpy from the point where the feed water enters a BWR to 
the point where the steam exits from the reactor consists of two stages. The water is 
first heated from its entering temperature to the T tem^ 
that is; tM saturation temperaiure> for the given system pressure. The water 
temperature does not rise above this value. The associated increase in enthalpy is 



This is the enthalpy of the saturated water. With the onset of boiling, the water 
absorbs an amount of heat equal to the heat of va/w/fea^^ 
unit mass of water that changes phase. Thus, the specific enthalpy of the steam is 
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= h in + c p (T)dT.rKft. 
Again, the rate at which heat is absorbedby the coolant is given by Eq. (8.4), 
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(8.5) 


where it? is the steam or feed water flow rate (the two are equal) in ib/hr or kg/hr. 

The advantage of expressing # in terms of enthalpy is that h has been com- 
puted using Eq. (8.3) or Eq. (8.5) and is available in tabular form for a number 
of substances. Abridged tables of h (these are called steam tables) are given in 
Appendix IV for water, helium, and sodium. It should be noted that only the dif- 
ference in enthalpy between two different temperatures is meaningful and useful; 
the absolute value of enthalpy has no meaning. Thus, the zero of enthalpy can be 
set arbitrarily, and the values of enthalpy computed from this starting point. In the 
case of water, for example, in the English system of units, h is taken to be zero for 
the saturated liquid at 32° F or in SI units as 0°C. 

Again, for water, it is observed in Appendix IV that the specific enthalpy 
is tabulated for saturated water and steam as a function of the temperature, and 
associated pressure in Table IV. 1, and as a function of pressure and the associated 
temperature in Table IV.2. (It must be remembered that there is only one pressure 
for a given temperature; conversely, there is one temperature for a given pressure at 
which water and steam are in equilibrium— i.e., at which water will boil.) However, 
the enthalpy of a liquid is not a sensitive function of pressure. Thus, the enthalpy 
of water in the pressurized water reactor is nearly equal to the value of h indicated 
in the tables for saturated liquid at the actual water temperature independent of its 
pressure. 2 This procedure is used in the first of the following two examples, which 
illustrate the use of 'Eq. (8.4) for PWR and BWR systems. 

Example 8.1 

A PWR operates at a thermal power of 3,025 MW. Water enters the reactor at 542.6°F 
and flows through the reactor at a rate of 136.3 x 10 6 Jb/hr. The system pressure is 
2,250 psia. 3 At what temperature does the coolant leave the reactor? 


2 Data for correcting saturation values of h for higher pressure are given in Table 4 of J. H. 
Keenan and F, G. Keyes, Thermodynamic Properties of Steam Including Data for . the Liquid and 
Solid Phases, New York: Wiley, 1959. In the temperature region from 400° F to 600 tt F and up to 
pressures of 300 psi, the correction for h amounts to less than J %. 

llSi^ns absolute-pressure in'lb/g£. The difference between the reactor pressure and at- 
mospheric pressure is called xauw pressure and is abbreviated psig. 
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Soluiion. According to Table 1.9, 1 k\V - 3,412 Bm/hr. Then 

q = 3,205 x tO 3 x 3,412 = 10, 321 x !0 6 Btu/hr. 
From Eq. (8.4), the rise in enthalpy through the core is 


- i - 10 ' 32i x l0 * 

?in ~ m'" 136.3 x I0 6 


75.7 Btu/lb. 


Interpolating from the values of enthalpy in Table IV. 1 for saturated water Jt J 
540° F and 550°F gives the value of h at 542.6°F as 539.7 Btu/lb. The enthalpy i>ff|§I 
emerging water is then 

h ml = 539.7 4- 75.7 = 615.4 Btu/lb. 

A second interpolation from Table IV 1 shows that this corresponds to water 1 
temperature of 599.1 °F. \Ans. \ 

Example8.2 . > _ c^V^? 

A BWR operating at pressure of>,6.895 kPa, produces 2.93 x 10 6 kg of stearoi 
hour. Feed water enters the reactor at0?0.6*<D. (a) What is the temperature of the 
steam? (b) At what power is the reactor operafitag? 

Solution. 

1. The steam emerges at essentially the temperature at which the water boilp! 
According to Table IV.2, at a pressure of 6.895 MPa, this isQf&8£?C. \Ans.] 

2, From Table IV.2, h of the steam is 2,773.2 kJ/kg. while from Table IV. I ^t 
1 90.6°C h of the feed water is 807.8 U/kg. Then Eq. (8.4) gives " : , • if 

<? = 2.93 x 1(^(2,773.2 - 807.8) - 5.758 x 10 6 J/hr = 1,599 MW. 


8.2 HEAT GENERATION IN REACTORS J 

The starting point in the design of a reactor cooling system is the deiermiriation 
of the spatial distribution of the heat produced within the reactor, lliis problenfe 
discussed in this section. ; f§| 

It was pointed out in Chapter 3 that the energy released in fission appears^ 
several forms — as kinetic energy of the fission neutrons, as prompt fission y-rays, 
as y-rays and $-rays from the decay of fission products, and in the emission of 
trinos. With the exception of the neutrinos, virtually all of this energy is uliimat|| 
absorbed somewhere in the reactor. However, because these various radiations ar. 
attenuated in different ways by matter, their energy tends to be deposited in djjgp] 
em locations. In the following discussion of the deposition of energy, it is assu " 


' neu- 
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(hk the reactor fuel is in the form of individual fuel rods. This is the case for virtu- 
ally all power reactor concepts except the molten salt breeder reactor (see Section 

lf\eat Production in Fuel Elements (Fuel Rods) 


As indicated in fable 3.6, the fissianjrag^ 168 ' 

ffvleV 'per 'fission! These highly charged particles have an extremely short range 
• (see Section 3.9), and therefore their energy is deposited near the site of the fis- 
sion within the fuel. Similarly, most of the 8 MeV of the fission product /J-rays is 
also deposited in the fuel. However, many of the /-rays froirtthe decaying fission 
products and those emitted directly in fission pass -out of the fuel since they are 
less strongly attenuated than charged particles. Some of these y-rays are absorbed 
;in the surrounding coolant and/or moderator, * n the thermal shield, or in the radi- 
ation shielding that surrounds the reactor. However, because of the proximity of 
fuel rods in most reactors, many of the y~rays am intercepted by and absorbed in 
neighboring rods. 

The prompt neutrons are emitted with; a total kinetic energy of about 5;MeV> 
■^l^ilfiflBnVln a thermal reactor, the bulk of this energy is deposited in the moder- 
ator as the neutrons slow down. The capture y-rays emitted following the absorp- 
tion of these neutrons in nonfission reactions are therefore produced and absorbed 
throughout the reactor. In a fast reactor, the fission neutrons slow down very little 
before they are absorbed, and their kinetic energy appears as an addition to the 
energy of the capture }/-rays. The delayed neutrons contribute negligibly to the en- 
^y^^gj^attor; |is already noted, none of the energy of the neutrinos is retained 
within a reactor. 

It should be clear from the foregoing remarks that the spatial deposition of 
fission energy depends on the details of the reactor's structure. NevertheiHss, for 
preliminary calculations, it may be assumed that approximately oneHhird of the 
total%-ray cner|ySaBout 5 MeV— is absorbed in the fuel. This, together with th^f 
168 MeV from the fission fragments and 8 MeV from "the /J-rays, gives 181 MeV' 
per fission (about 90% of the recoverable fission energy), which is deposited in^ 
thesftiel, ; m^ v '£MJ^9| l^£ ^ ssicltn site. The remainder (about" 

20iMeVf of the recoverable energ^ trv 

various structural materials, and in the blanket, -reflector, and shield. * 

- 

The rate at which fission occurs in the fuel, and hence the rate of production 
of heat, varies from fuel rod to fuel rod, and it is also a function of position within 
any given rod. In particular, if £j is the energy deposited locally in the fuel per 
fis$jort then the rate of heat production per unit volume at the point r is given by 
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where S /r (E) is the macroscopic fission cross-section of the fuel and $|r, |ff 
the energy-dependent flux as a function of position. T^d r naturaJ units of |JF 
|^^fe^.H5m 3 ; However, for engineering calculations, q' n should be converted tcj . 
unitsof Btu/hr-ft 3 or kW/lit using the conversion factors in Appendix 1 c t 

In the thermal reactor, most of the fissions are induced by thermal neutpSt 
In this case, according to the results of Section 5.9, the integral in Eq. (8.6) 
written as 


where l/fcisteth^ 

Mx% In the framework of a multigroup calculation, Eq. (8.6) can also be express! 
in the form 

where the sum is carried out over all the groups and the notation is the same as in 
Section 5.8. : | 

The spatial dependence of the flux depends on the geometry and struc^ 
of the reactor, However, many heat removal calculations are carried out for the 
theoretical case of the bare, finite cylinder. For a thermal reactor, the thermal ffixjjt 
is then (see Table 6.2) 

. .. ^ 3.63 P , L2M5r\^Jmz\ , ^fJ<a*s 

E R T f v J ? \ h J ■ y " W 

where P is the total power of the reactor in joules, E R is the recoverable energ)Jj|r 
per fission in joules, V is the reactor volume in cm 3 , and R and H are its bif w 
dimensions in centimeters to the extrapolated boundaries. In obtaining Eq..(8-' 
it was assumed that the fuel is homogeneously distributed throughout the reactor| 

4 The following notation for heat production and flow is well established in heat transfer 
ature and is used in this chapter: 

power density, rate of energy production per unit volume; BTU/hr-ft 3 or kW/HU 
heat flux; BTU/hr-ft 2 or W/cm\ 
linear power: Btu/hr-ft or kW/m, 
heat transfer rate; Btu/hx or MW. 


r 



hr~ft 3 or kW/lit, *\ 
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and L / is the value of the macroscopic fission cross-section of this mixture in units 
ofcm~ ! . 

Equation (8.9) can also be used to approximate the fluxjn a reactor where the 
fuel is contained in separate fuel rods, provided the value of £ / is computed for the 
equivalent homogeneous mixture. For instance, suppose that there are n fuel rods 
of radius a and length //, the height of the core. Then if £ /r is the macroscopic 
fission cross-section of the rod, the total fission cross-section in the entire core 5 is 
E/> x nira l H. Therefore, the average value of E/ in the core is 


t - ^frnJta r H _ E fr ncr 


(8.10) 


and the flux is 


_ 3.63P/? 2 f /2.405A / nz\ 

E H E fr Va 2 n \ R J \H J 



(8.11) 


where V = n R 2 H has been substituted. 

When the expression for 4> T in Eq. (8.1 1) is introduced into Eq. (8.7), the rate 
of heat production per unit volume of a fuel rod becomes 


u , r Jn cos ~ ) . 


(8.12) 


It should be noted that, in this procedure for obtaining q m , any variation of the flux 
across the diameter of the fuel rods has been ignored. The dependence of q m on r 
in Eq. (8.12) gives the change in the flux from rod to rod across the diameter of the 
core — not across any individual fuel rod. Indeed, given the derivation of Eq. (8.12), 
q m must be taken to be constant across individual rods. This does not introduce 
significant errors in heat transfer calculations, especially for the small-diameter, 
weakly absorbing (low-enrichment) fuel rods currently used in most power reac- 
tors. 

From Eq. (8. 12), it is evident that the maximum rate of heat production occurs 
in the middle (z = 0) of the central rod (R = 0). In this case, both functions in Eq. 
(8,12) are unity, and the maximum value of q m is 

(8.13) 


M6PE, 



5 Recall chat macroscopic cross-section is cross-section per unit volume. 
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TtTe maximum me of heat produc r 0 is 


/" <» / 2:405r \ 


(8.14) 


The total rate at which heat is produced in any fuel rod is given by the integral I 


q r (r) ss ?ra 


up. 

2 1 q m {r<z)dz. 


fN/2 

J-H/2 


Introducing q"'{r, z) from Eq. (8.12) gives 
lA6nPE d 


_ 23TPEJ • Y 2^05r \ 

. .. .' jr ■."V / 


(8* I'll 


The formulas derived previously for the hare cylindrical reactor should not 
be taken too seriously for calculations of the heat production in a real reactor, In v A 
particular, Eq. (8.13) considerably overestimates the value of for a reflectgd^ ^j 
and/or nonuniform!}' fueled reactor, which has a smaller maximum-to-average flux^ ffg|| 
ratio than a bare reactor (see Section 6.3). To see this effect, it is first noted that / 

where 0 max is the maximum value of the thermal flux. The total reactor power is ^ 
given by £ ^ ♦ 


5 ^ \ 


\ 


where <f> w is the average thermal flux and E/ is again the macroscopic fissici||!|^ 
cross-section averaged over the entire core volume V, Dividing Eq. (8.1.6) by'i||pM 
(8. 1 7) and rearranging gives | | ] 

q ™~ Z f E R <p m V ~ E f E K V ■ mi ] 

where Q is the maximum-to-average flux ratio. Finally, introducing Eq, (8.10) lor 
1/ yields ~ i 


m -m^H^ 

£ ZWM*^.. 

m wm:-- 


9max ~ ahiVEk ~~ 


(8,19) >: 

Suppose that Q = 2.4Pa?reasonable value for an actual reactor. Then comparing 1 ™ ■ 
Eqs. (8. 1 3) and (8.19) shows that " 


' ■ 1 
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Example 8 J 


The extrapolated dimensions of a certain pressurized water reactor are = 67 in and 
H = I44in. The reactor operates at the thermal power of 1 ,893 MW. It contains 193 
fuel assemblies, each consisting of 204 UO2 fuel rods 0.42 in. in diameter. Assuming 
that the assemblies are uniformly distributed throughout the reactor, calculate the 
total energy production rate and the maximum energy production rate per ft 3 in rods 
located (a) at the axis of the reactor, (b) 20 in from the axis. 

Solution. 

1. There are n = 193 x 204 = 39372 rods in the reactor. Assuming that E (J = 
1 80 Me V and E R = 200 Me V, Eq. (8. 15) gives for the rod at r = 0 

, m 2.32 x 1893 x 180 A <#mMf 
* (0> -- 39.372 x 200 " ai00MW 

= 3.43 x 10 5 Btu/hr. [ito.] 
Comparing Eqs. (8. 13), (8. 1.4), and (8.15) shows that 


2 Ha 


Therefore, 


3.43 x 10 5 


= 4.66 x 10 7 Btu/hr-ft 3 (482.3 kW/lit). [Ans.\ 

2. To obtain similar results for a rod located at r = 20 in, it is merely necessary 
to multiply the prior answers by the factor appearing in Eqs. (8.14) and (8.15), 
yielding, ■> 


Jo 


( 2.405/- \ /5L405 x 20\ 


Thus, for this rod, 


q r (20 in) = 3.43 x 10 s x 0.875 

= 3.00 x )0 5 Btu/hr (87.9 kW) [Am.) 
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m 

^(20 in) = 4.66 x 10 7 x 0.875 || 
s 4.08 x 1 0 5 Btu/hr-ft?(4.2 kW/lit). [Arts.] } \ 


Radiation Heating 

As noted already, roughly 10% of the recoverable energy of fission is absorbed 
outside the fuel. In thermal reactors, the kinetic energy of the fission neutrons is 
deposited in the surrounding moderator and coolant in more or less the same spati?* 1 - 
distribution as the fissions from which these neutrons originate. However, only 2| 
to 3% of the fission energy appears in this form, and it is often assumed that thii 
kinetic energy is deposited uniformly throughout the core. 

The calculation of the energy deposition from the longer range x-rays .is 
a more difficult problem. In principle, this can be determined by evaluating the, 
integral 


q'"{r) - j <j>y(r.E y )Eyii a (Ey)dE y , 


(8.20)' 


where <p y (r T E y ) is the y-ray flux as a function of position and energy, and n a {E Y ) 
is the linear absorption coefficient. The computation <f> y (r,E y ) is complicated by 
the fact that y-rays undergo multiple Compton scattering— a problem discussed in 
Chapter 10. However, in source-free regions outside the core, such as the thermal 
shield or reactor vessel the y-ray flux may be presumed to fail off approximatejS 
exponentially, provided the region in question is not too many mean-free patH|| 
thick. In this case, 


q m ir) = f -e~»* iE r )r <f>yo(E y )E Y ix a (E r )dE y . 


Here, (f> y o(E Y ) is the y-ray flux as a function of energy incident on the regioii^ 
Equation (8.21) is usually evaluated by dividing the y-ray spectrum into a numbe^f 
of energy groups. This procedure is illustrated in Section 1 0. 1 1 . f ;% 


energy groups, I his proced 
Fission Product Decay Heating 


After a few days of reactor operation, the ft- and y-radiation emitted from decaying: 
fission products amounts to about 7% of the total thermal power output of tfe 
reactor. When the reactor is shut down, the accumulated fission products continu| 
to decay and release energy within the reactor. This fission product decay energy 
can be quite sizable in absolute terms, and a means for cooling the reactor : cO.?i 
after shutdown must be provided in all reactors except those operating at very \om 
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trons in the m U. must also be taken into account. Using the equations of radioac- 
tive decay derived in Section 2.9 ? it is easy to show that 


-~ = 2.28x l(T 3 C 


e ~4.9l x 10~"0 , — 4.9I x lO -4 * 


and 


"39 


= 2.I7 x 10 


\°nsJ 


id 


. e -iA\ x It)"** ^-3.4 1 x lir''. 



-7.0 x I0~ 3 (1 -e- 4 -» lxl0 " 4 ' 0 )^ ,x '0" 4 "]. 


(8.24); 


In these equations, P20 and P 39 are the decay powers of 239 U and 239 Np, respec- 
tively. C is the conversion factor for the reactor, and W a2 s and a f2 s are the effecttyf 
thermal cross-sections of 235 U. The times / 0 and t s are again in seconds. 

Example 8.4 j, ;sjfj 

A certain -"U-fueled reactor operates at a thermal power of 825 MW for 1.5 years 
and is (hen shut down, (a) Using Fig. 8.3, compute the decay energy at the following 
times: at shutdown, 1 hour after shutdown. 1 year after shutdown, (b) If the conver- 
sion factor for the reactor is C = 0.88, what are the contributions to the decay energg 
at the prior times due to and 239 Np? 

Solution. . ( fli 

1. In this problem, P {> = 825 MW and % = 1 .5 years = 1 .5 x 3. 16 x 10 7 secj 
Then at shutdown, t s =10"' sec (the shortest time in Fig. 8.3). and from £j§ 
(8.22). 

P P(W~ l ) F(4.74 x IP 7 ) 


Po 


Po 


From Fig. 8.3, this is 


— = 0.070 - 0.0007 ~ 0.070, 
"0 


and P = 825 x 0.070 = 57.7 MW. \Ans.\ For t, = I hr = 3.6 x 10 3 sec: 

L - P(36 * l0l) _ p (4.74 x 10 7 ) 
Po ~ P 0 To 

= 0.014 - 0.0007 ~ 0.014. 

so that P = 825 x 0.014 = 11.5 MW. [Arts.] Finallv. for i t = I vr = 3.16 % 
10 7 sec. f 0 + i, =7.90 x I0 7 sec: 
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p _ P0M * l0?) - P(1 - 90 * l0?) 
To" Pa Pa 

= 0.00079 - 0.00063 = 0.00016. 
P = 825 x 0.00016 = 0.132 MW. [Arts.] 

* 

2. if the difference between the non-l/u factors for a tt2 5 and 0725 is ignored. 
Van/ops = 681/582. Then with / 0 = 4.74 x 10 7 sec, the exponential in the 
brackets of Eq. (8.23) is negligible, and with t s = 0 this equation gives 

^ = 2.28 x 10- 3 x 0.88 x $|I = 135 x I0~ 3 . 
Pi) 582 

• Then /fc - 825 x 2.35 x 10~ 3 = 1 .95 MW. f/W.J 
For t $ - 1 hr, 

% s 1.95 exp[~4.9 x !0~ 4 x 3.600] = 033 MW. [Ans.\ 

With t s = I year. / J 29 ~ 0, as the half-life of 239 U' is only 23.5 minutes. [4/i.v.J 
For P3<>, using Eq. (8.24) in a straightforward way. the results are: at shutdown, 
j .84 MW; after I hour, ^ 1 .84 MW; after I year, ~ zero. \Ans.) 


3 HEAT FLOW BY CONDUCTION 


Energy is removed from a reactor by two fundamentally different heat transfer 
processes — conduction and convection. In conduction, heat is transmitted from one 
location in a body to another as a result of a temperature difference existing in the 
body — there is no macroscopic movement of any portion of the body. It is by this 
mechanism, as shown in this section, that heat produced in a fuel rod is transferred 
to the surface of the rod. Heat convection involves the transfer of heat to a moving 
liquid or gas, again as the result of a temperature difference and the later rejection 
of this heat at another location. Thus, the hear conducted to the surface of a fuel 
rod is carried into the coolant and out of the system by convection. Such convective 
heat transfer is discussed in Section 8.4. 

For completeness, it should be mentioned that heat can also be transferred as 
thermal radiation across a vacuum or other rarefied space between a hotter body 
and a colder one. However, this process is of relatively little importance except in 
some gas-cooled reactors and is not considered further. 

The Equations of Heat Conduction 

The fundamental relationship governing heat conduction is Fourier's taw, which 
for an isotropic medium is written as 

q* = ~k grad T. (8.25) 
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Here q", which is called the heat flux, is defined so that q" • n is equal Jfil| 
rale of heal flow across a unit area with unit outward normal n. Thus, q" is enhiefe^i' 
analogous to the neutron current density vector defined in Chapter 5. In the Bridal 
system, q" has units of Btu/hr-ft 2 and W/m 2 in SI units. The parameter k : in fL 
(8.25) is called the thermal conductivity and has units of Btu/hr-ft-°F in the English 
system and W/m-K in SI units. Values of & for a number of important suosttllfl 
are given in Appendix IV. In general, k is a function of temperature. The function' * 
T in Eq. (8.25) is the temperature in °F or K as appropriate. The similarity between 
Eq. (8.25) and Fick's law of diffusion should be especially noted. I 
Consider an arbitrary volume V of material throughout a portion of which 
heat is being produced. From the conservation of energy, the net rate at which heat 
flows out of the surface of V in the steady state, must be equal to the total rale at?,,,, 
which heat is produced within V. If this were not the case, the substance wouU|j^| 
change temperature and therefore would not be in a steady state. In equation forraF* 
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term appearir 
as neutrons ai 
When F 
equation is di 


, This is called 
I Poisson'seqi 


(8.3 1) reduc* 


Net rate of flow 
of heat out of V 


rate of heat production 
ithln V 


>]=„. 


The net rate at which heat flows out of the surface of V is 


Heat flow 


/ q" n 

J A 


dA, 


(8.2il 


(8.27) 


IIS 


where the vector q" is as defined earlier, n is a unit vector normal to the surtpE 
and the integral is taken over the entire surface. From the divergence theorem, j 
(8.27) can also be written as r 4 

Heat flow = / div q'V/V, 

Jv • 3Sj 

The total rate of heat production within V is equal to 
Heat production = / q"'dV, 


where q'" is the rate at which heat is produced per unit volume. % ;| 

Equations (8.28) and (8.29) can now be introduced into Eq. (8.26). Since the ^ 
integrals are over the same arbitrarily selected volumes, their integrands mu$t$f 
equal and the following expression is obtained: f 

divq*-f w = 0.. ' 1 


which is call 
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This result is the steady-state equation of conductivity for heat transfer and is|pg 
gous to the equation of continuity (Eq. 5.15) discussed in Section 5.3 in conn|c|^ 
with neutron diffusion. There is no term in Eq. (8.30) equivalent to the absorpjg 
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term appearing in the neutron case; heat simply does not vanish within a medium 
as neutrons are inclined to do. 

When Fourier's law, Eq. (8.25), is substituted into Eq. (8.30) and the resulting 
equation is divided through by k. which is assumed to be constant; the result is 


v 2 r= V = 0. 


(8.31) 


This is called the steady-state heat conduction equation and is of a form known as 
Poisson's equation. In a region where there are no heat sources. q m = 0 and Eq. 
(8.31) reduces to 


V 2 T = (\ 


(8.32) 


which is called Laplace's equation. 

These results are now applied to some problems of interest in nuclear reactors. 
One of the central problems, as is seen, is the calculation of the heat that can be 
transferred out of a fuel rod and ultimately into a coolant for a given maximum 
temperature in the fuel. The maximum fuel temperature is a preset condition that 
must not be exceeded for reasons of safety. 

Plate-Type Fuel Elements 

Consider first a plate-type fuel element or rod like that shown in Fig. 8.4, consisting 
of a fueled central strip (the "meat") of thickness 2a surrounded on all sides by 
cladding of thickness b. It is assumed that heat is generated uniformly within the 
fuel at the rate of q m Btu/hr-ft 3 or kW/lit, and that the temperature has reached a 
steady-state distribution throughout the element. 

Ordinarily, the total thickness of such an element is small compared with 
either its width or length. It is reasonable to ignore the negligible amount of heat 
flowing out through the edges or ends of the element. In short, the heat flows only in 
the x direction, where x is the distance from the center of the element normal to the 
surface (see Fig. 8.4). The temperature distribution in the fuel is then determined 



Figure 8.4 A plate-type fuel 
element. 
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by Poisson 's equation in one dimension, 

d 2 T q>" I 

7^ + IJ =0 ' (8-33 

where *y is the thermal conductivity of the fuel. 1 
Two boundary conditions must be specified with a second-order different 

equation like Eq. (8.33). In the present problem, these are M 


(i) T(0) = T a; 

where T,„ is the maximum (central temperature) in the fuel, and 

(IT 
(ii) — = 0 

ax .Hj 

at v = 0. The second condition follows from the symmetry of the problem which 
permits no heat flow at the center of the fuel. ' 

Integrating Eq. (8.33) twice gives the general solution . j 

_ cf , 

where C, and C 2 are constants. Placing a- = 0 immediately gives C, = T,„ Als< 
in view of condition (ii), C, must be taken to be zero. The temperature within 
fuel is therefore 


T — T 


r , 
57'" 


Using this equation, the temperature T, at the surface of the fuel (at the. 
cladding interface) may be evaluated. Thus, writing x = a in Eq. (8.34) give! 


T S = T„ 


q"'a 2 
2k f ' 


1 he total rate of heat production in the fuel is equal to q'" multiplied ^ 
fuel volume 2Aa, where A is the area of one face of the fuel. In die steady state! 
of the heat produced wi thin the fuel flows out of the fuel. The heat flowing i 
one face of the fuel is therefore 


q = q"'Aa. 


This result can also be obtained from Eq. (8.34) using Fourier's law. Thus, 
flux (heat flow per unit area) is given by 


if . 

q = k f 


sir 
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where the derivative is to be evaluated at x = a. Inserting Eq. (8.34) and perform- 
ing the differentiation yields 

q" = q"'a. 

The total rate of heat flowing out of one side of the fuel is then 

q = q"A = q"'Aa, 

which is the same as Eq. (8.36). . ,„ . 

It is sometimes convenient to rewrite Eq. (8.36) ,„ a form y*£^J" 
been eliminated through Eq. (8.35). Thus, solv.ng Eq. (8.35) for q and .nsertmg 
the result into Eq. (8.36) gives 

= T "' ~ r * (8-37) 
q a/2k/A' 

This expression may be viewed as the heat transfer analogue of Ohm's law in 
electricity— namely, 

(8.38) 


V 


where / is the current V is the potential difference, and R is the resistance In the 
;^t casf Responds to ' T m - % is analogous to the potential dvfterenc. 
ZTalKfA I called the thermal resistance. The value of (hi. analogy » clear 

Sh ° rt Turning next to the temperature distribution in the cladding of the plate-type 
fuel elemenfshown in E.g. 8 4, it is first observed that, since there ts htt e or no 
heat genm"ed in this region, q»> = 0 and the heat conduction equauon reduces to 

^1 = 0 (839) 
dx 2 

This is the simplest form of Laplace's equation (Eq. 8.32). The boundary conditions 
are now 

(i) T(a) = T s 

(ii) T{a + b) = T c , 

where T c is the temperature at the outer surface of the cladding. 
Integrating Eq. (8.39) twice gives 
*" T = C\x+C 2 . 

i The constants C, and C 2 are easily found from conditions (i) and (ii). and the final 
expression for T is 
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or 



(8.42) 


This formula gives the rate at which heat flows through one side of the fuel plate in 
terms of the difference in temperature between the center and surface of the plate. 
In the electrical analogue, T m - T c is the difference in potential and 

R = *h + A » 43 > 

2k f A k c A 

is the total thermal resistance of the fuel and cladding. 

Equation (8.43) shows that, when heat passes through a succession of two 
materials, the situation is analogous to an electrical potential across two resistances 
connected in series, where the total resistance is the sum of the two. This conclusion 
is readily generalized to a sequence of any number of materials — the total thermal 
resistance is the sum of the resistances of each material 

In the derivation of Eq. (8.42), it is assumed implicitly that the fuel and 
cladding are tightly bonded together at their point of contact, Some fuel elements, 
however, have a thin region of bonding material between the fuel and cladding, 
whereas others have a narrow region of gas between the two. The total thermal 
resistance is larger in this case and is obtained by adding an appropriate resistance 
term to Eq. (8.43). The temperature difference T m - T c must then be greater for a 
given heat flow out of the surface of the fuel element. 

Cylindrical Fuel Rod 

Consider a long cylindrical fuel rod of radius a surrounded by cladding of thickness 
b. It is assumed again that heat is produced at the constant rate q m within the rod 
and that there is no heat released in the cladding. 

The temperature in the rod is only a function of the distance r from the axis 
of the rod, so that in cylindrical coordinates the heat conduction equation is 

d 2 T IdT q m n , OAAS 
77 + -T + T-=0. (8.44) 
dr 2 r dr kj 

The boundary conditions appropriate to this problem are 

(i) T is nonsingular within the rod, 
(u) T(0) = T mt 

where T m is the central temperature of the fuel. 
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ll can be readily verified by substitution that the general solution to Eq. ( Jf 


„< r i 

r = -I-_ + c, Inr + CS. 

where C, and C 2 are constants to be determined. In view of the boundary condi- - ■ 
lions, it is evident that C, = 0 and C 2 = T m . Thus, the temperature within the roM 
is 


T = T — 


q"'r 2 
4k f ' 


. 11 
(8.45) . r 


The rate at which heat is produced within the rod, and therefore the rate at ' 
which it flows out of the rod, is equal to 

q=na 2 Hq'", 

where H is the length of the rod (approximately equal to the height of the core). 
Solving for*/"' and introducing this into Eq. (8.45). evaluated at /• = a where the M 
temperature is T $ , then gives 



q = 


\/(4xHk f y 

In view of Eq. (8.47), the thermal resistance of the fuel is clearly 

»,-- U 

1 4irHkj 
The heat conduction equation for the cladding is 


d 2 T 1 dT 
_. + -__ = 0 . 
dr- r dr 


which has the solution 


T = C, In r + C 2 . 


Using the boundary conditions 

(i) T(a) = T s , 

(ii) T{a + b) = T c 

to determine the constants C, and C 2 gives 

?; ln(r/ + />) - T c In a - (T s - f,) In r 


T = 


ln(l + 
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From Fourier's law, the total heat flowing out of the cladding is 
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q = -In (a +b)Hk c 


d7' 


; 


||%valuated at r = 0 + Carrying out the differentiation of Eq. (8.49) yields 

2xHkc(T,-T e ) 


q = 


, *!/The thermal resistance of the cladding is therefore 


ln(1 +b/a) 
I is therefore 
= ln(l +f>/o) 

#1:'.. 2/rWAv " 

• Frequently, b is much less than a. In this case, since 


(8.50) 


(8.51) 


In 


KB 


III" Eq, (8.5 1 ) can be written as 


1M 


R c ^ 


2naHk c 


(8.52) 


In terms of the overall temperature difference between the center of the fuel 
and the surface of the cladding, the heat flowing out of the rod is 


Rr + R</ 


(8,53) 


where R f is given by Eq. (8.48) and R c is given by Eq. (8.5!) or (8.52). This result 
may be verified by solving Eq. (8.47) for T m - T s and Eq. (8.50) for T s - T c and 
adding the two expressions. 

Space-Dependent Heat Sources 

In both of the foregoing calculations of the temperature distribution and heat flow 
in plate-type and cylindrical fuel rods, it is assumed that the heat is produced at 
the constant rate q m throughout the fueled portion of the rod. Accordingly, the 
temperature doe's not vary along the length of the rod, and it is possible to relate 
the total heat flow out of the entire rod to the difference in temperature between the 
center and surface of the rod. 

However, it is shown in Section 8.2 that q m does, in fact, vary approximately 
as (see Eqs. [8.12] and [8 13]) 
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where z is measured from the midpoint of the rod. There is negligible variation 
in q'" in the x and y directions of the plate-type element or in the radial direction 
of a cylindrical rod. A nonuniform heat distribution like that given in Eq. (8.54) ; 
would be expected to give rise to a nonuniform temperature distribution along the 
length of a fuel rod— and. indeed, to some extent this is the case. In a power reac- 
tor, however, the coolant passes along the fuel rods at a temperature considerably 1 
lower than the temperature of the fuel, and the temperature gradient is therefore 1 
much steeper across the diameter of the fuel than along the length of the rod. As a 
consequence, most of the heat generated in the fuel flows directly to the surface of 
the rod; there is little or no heat flow in the z direction. 

It is possible to generalize the calculations made in this section to the case 
where q'" is a function of z. This can be done by dividing the calculated value of |, 
which is the the total heat flow out of the surface of a fuel rod, by the surface^ area | 
of the rod. This gives the heat flux, q", which is independent of position if q'" is 1 
constant, but varies with position if q'" is not. For instance, if Eq. (8.46) is divided 3 
by 2n(a + b)H, the area of a clad fuel rod, the result is 


2(a + b) 


(8.55) 


This is the heat flux at the surface of the rod as a function of z. Similarly, q"{z) 
can he found as a function of the temperature difference T m - T c , where' bc-ufl 
temperatures may be functions of 2, by dividing Eq. (8.53) by 2n(a + b)H: 


T m (z)-T c (z) 
27t(a + b)H{R f + R t ) 


(8.56) 


Example 8 S 

The fuel rods for the reactor described in Example 8.3 consist of a fueled ponm 
0,42 in. in diameter that is clad with Zircaloy-4, 0.024 in thick. Each rod is 1 2 ft lonj 

Given that the center temperature of the fuel is 3.970 °F at the midpoint of the centra. 

rod, calculate at this point the (a) heat flux out of the rod; (b) outer temperature of:; 
the cladding. 

Solution. 

I. From Example 8.3, q m at the midpoint of the centra] rod is 4.66 x 10 7 BTU/h^; ;| 
ft 3 . The radius of the fuel is 0.21/12 ^ 0.0175 ft, whereas the cladding thick' 
ness is 0.024/ 1 2 = 0.002 ft. Then using Eq. (8.55), 
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(0.0175) 2 
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x 4.66 x 10' 


2(0.0175 + 0.002) 
= 3.66 x 10 5 Btu/hr-ft 3 
= 1 15.5W/cm 2 . \Ans.\ 

2. From Table IV.6, k f = 1.1 Btu/hr-ft-°F and k c ~ 10 Btu/hr-ft-°F. Equation 
(8.48) then gives for the thermal resistance of the fuel 


I 


4ir x 12 x 1.1 
For the cladding, using Eq, (8.51), 


6.03 x J0" J *F-hr/Btu. 


ln( 1 +0.024/0.21) 

R c = — L ' = 1.43 x 10~ 4 T«hr/Btu. 

2?r x 12 x 10 

The total resistance is then 6. 17 x 10 3 °F~hr/Btu. 

Introducing the prior parameters into Eq. (8.56) then gives 

T c = 3,970 - 2^(0.0175 + 0.002) x 12 x 6 17 
x Hr 3 x 3.66 x 10 5 = 650°E [yln*.] 

Exponential Heat Sources 

It is often necessary to calculate the temperature distribution and heat transmission 
m reactor shields and pressure vessels in which radiation energy is deposited more 
or less exponentially. Consider a slab of thickness a whose surfaces are held at the 
temperatures T x and T 2 , as shown in Fig. 8.6. If x is measured from the surface as 



/-nx) 











s 




■1 


Figure 8,6 Slab containing 
exponentially distributed heat 
source. 
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indicated in the figure, the heat source distribution is given by 


(8|| 


where S and \x are constants. For y-rays, /< may often be taken to be the lineal 
absorption coefficient. 

Introducing Eq. (8.57) into the one-dimensional heat conduction equatio 
(Eq. 8.33) gives 

d2 ' r s „. „ 

The general solution to this equation is easily found to be ! 


kix 2 


i':H 


where the constants C } and C 2 are found from the boundary conditions 

r(0) = r l? 7» = r 2 . 

The final expression for the temperature distribution is 

This equation is the starting point for calculations of heat removal frol B'i^^L| 
which radiation is absorbed exponentially. It is also used to compute thie^aMt 
stresses that accompany such nonuniform temperature distributions. As is sh 
in the problems at tJhe end of the chapter, under certain conditions the temperai 
in the slab may rise to a maximum value that is larger than either T\ or T 2 . 
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T = Tt + (T 2 - r,)- + ~ \l - e~*" - -<1 - • ' M 

a kii 2 L a J 
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8.4 HEAT TRANSFER TO COOLANTS 

As noted earlier, the heat produced in the fuel or deposited by radiation : i| 
parts of a reactor is transferred to a coolant of one type or another, which, ir, ~ 
carries the heat outside the system. The fundamental relation describing the:tn||| 
of heat from a heated solid to a moving fluid (liquid or gas) is Afew%?^^ 
cooling: 

In this expression, q t! is the heat flux in Btu/hr-ft 2 -°E T c is the temperatii^^ffl 
surface of the solid (in the case of reactor fuel this is the outer temperature of^al 
cladding), and T h is an appropriate reference temperature of the fluid. 
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The numerical value of the heat transfer coefficient depends on many factors, 
including the nature of the coolant, the manner in which it flows by the heated 
surface, and the coolant temperature. Methods for calculating h from experimental 
data are presented later in this chapter. At this point, however, it is of interest to 
note the range of values of h that are applicable to the coolants in power reactors. 
For ordinary and heavy water, h usually lies between 5,000 and 8,000 Btu/hr~ft 2 -°F 
(25 and 45 kW/m 2 K); for gases, h ranges between about 10 and 100 Btu/hr-ft 2 -°F 
(55 and 550 VV/nr K); and for liquid sodium, h is generally between 4.000 and 
50,000 Bm/hr-ft 2 ~°F (20 and 300 kW/m 2 K). 

In most reactor heat transfer problems, the fluid flows along the fuel in well- 
defined coolant channels. Since the temperature of the fluid ordinarily varies with 
position across each channel it is possible to define the temperature T h in Eq. 
(8.60) in any number of ways. In most heat transfer calculations, T b is taken to be 
the mixed mean or bulk temperature of the fluid, This is defined by the formula 


fpc»vTdA e 

Tt= J t P , (8.61) 

}pc p vdA c 

where p is the fluid density, c p is its specific heat, v is the fluid velocity, T is the 
temperature, dA c is the differential of channel area, and the integrations are carried 
out over the cross-section of the channel In general, each of the variables in the 
integrals may be a function of position across the channel The temperature T b 
defined by Eq. (8.61.) is the temperature the fluid would achieve if it were allowed 
to mix adiabatically — for instance, if the fluid flowing out of the channel were 

icollected and mixed in an insulated container. Measured values of ft are normally 

-given in terms of the bulk temperature of the fluid; 

i„ The total rate of heat flow across an area A between a solid and a fluid is 

q=q"A^hA(T c -T h ), (8.62) 
. Written in the form of Ohm's law, Eq. (8.62) is 

: q = -777^* (8.63) 

if ha 

The denominator in this equation is the thermal resistance for convective heat trans- 
fer, 

ft-jL. (8.64) 

It is now possible to return to the fuel element problems considered in the 
last section and compute the rate at which heat is transferred to a coolant, for a 
given difference in temperature between the center of the fuel and the fluid. For 
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the cladded plate-type fuel element, the total thermal resistance is the sum of f ' 
(8.43) and (8.64): ,M 

.. lilll III 

a b 1 ■^pMi.iip! 

2A V ^ + M + M' : ./ '* 

and the total heat flow through one side of the fuel element to the fluid is 


llll 


T m — Th 


q = 


a b 1 

2F^4 + M *~hk~ 


Similarly, from Eqs. (8.48) and (8.51). the total thermal resistance for a cladcyl&jj 
drical fuel rod is 

R _ I ln(l + b/a) J_ 
~ 4izHk f + 2nHk c + hA' 


where A 


2n(a +b)H. Kb 

I 


R = 


4- 


+ 


I 


AnHkj IrtallK. hA 
In either case, the rate of heat flow into the coolant is 

T m — Tb 


R 

These results, which apply to the case of uniform heat production along the; 
length of the fuel, can be generalized to the more realistic situation where q" ^f||jj 
pends on z. Thus, as explained in the preceding section, most of the heat produced 
in the fuel flows directly to the coolant in a direction normal to the axis of the W { - 
rod. The expression for the heat (lux. obtained by dividing q in Eq. (8.69) by the 
surface area of the rod, is therefore valid at every point along the rod. Thus 
cylindrical rod, it follows that 



q"(z) 


T m (z) - T„(z) 


2n(a + b)HR 
In terms of T c . the outer temperature of the cladding, 

T c (z) - r b (z) 


q"(z) = 


27T(a+b)HR h 
= h[TAz) - T„(z)h 


which is simply Newton's law, Eq. (8.60). 
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The fuel rods in Examples 83 and 8,5 are cooled with pressurized water. Given that 
the heat transfer coefficient is 7.500 Btu/hr~ft 2 ~°F, calculate the bulk temperature of 
the water opposite the midpoint of the hottest rod. 

Solution According to Example 8.5, at the point in question, q' 1 = 3.66 x 10 5 
Btu/hr-ft 2 and T c s 650°E From Eq. (8.60). 

T b = 650 ~ 3,66 x 10 5 /7.500 

: - =60l°F.tAw.] 

■::| For reasons discussed later, this is not the highest temperature of the coolant al- 
| though q m is greatest at this point. 

fTemperature along a Coolant Channel 

JA$ the coolant moves along the fuel it absorbs heat; as a result, its temperature 
continually increases, However, the temperature does not increase at a constant 
frate since the heat is released from the fuel nonuniformly — according to 


9 =% 


cos 


")■ 

HI 


(8.71) 


This equation applies to the central (r = 0) rod; for other rods, q m must be mul- 
tiplied by J(}(2AQ5r/R). A simple method is now discussed for estimating the 
variation in the temperature of the coolant and the fuel as a function of z, when the 
heat is produced in the fuel according to Eq. (8.71). 

Before proceeding, it should be noted that, except for the fuel rods at the edge 
of the core, there is only one coolant channel associated with each fuel rod. This 
may be taken to be either the volume between neighboring fuel rods, with a portion 
of each rod contributing heat to the channel, or the volume of coolant in a unit cell 
surrounding one rod. As shown in Fig. 8,7 for a square array of rods, the volumes 

tjbf coolant associated with the surface area of one rod are the same in both cases. 
Consider a coolant slab of thickness dz as it moves along a coolant channel. 
The volume of this slab is equal to A t dz, where A c is the cross-sectional area of 

I'the channel, and its mass is pA c dz* where p is the coolant density. As the slab 

F progresses up the channel the distance dz, it absorbs heat from the fuel, which 
raises the temperature pf the slab the amount dT*. 6 The heat required to produce 

|the change dT b is equal to 

;| pA c dzc p dT b , 



i 6 Throughout this discussion, the coolant is assumed not to boil. If it is boiling, then the ab~ 
'■' sorption of heat goes into vaporization of the coolant at constant temperature, 
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Figure 8.7 Equivalent cooling 
channels in a square fuel lattice. 


where c p is the specific heat per unit mass. The rate at which heat is added ®! 
flowing coolant is obtained by dividing this quantity by the time dt required fp| 
coolant to move the distance dz. Thus, 

dq pA^CpdTt, = pA c vc p dT h , 

where v is the velocity of the fluid. However, pA - cv is equal to the rate of § 
of the coolant through the channel, w, and so 



dq rr WC p dTh 


As noted earlier, most of the heat produced in the fuel flows directly into the 
coolant and not along the interior of the fuel Therefore, the heat dq in Eq. (8.72] 
is generated in a fuel volume of A f dz. where A f is the cross-sectional area of the 
fuel portion of a fuel rod. Thus, 

WCpdTh = q m 'Afdz. 

Introducing the function q m from Eq. (8.7 1 ) and integrating from the point of entry 
of the coolant at z = ~~H/2, where its temperature is T Hh to an arbitrary point : 
along the channel gives 


(8.72) 



7* , 7" * 'mar I 


7i we 


I + sin 


JlWCn 


P I 
I + 


-mi 


where V f is the volume of the fueled portion of the rod. 
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Equation (8.73) gives the temperature of the coolant as a function of position 
along the central or hottest channel The temperature along other channels can J>e 
obtained by multiplying the second term in Eq. (8.73) by the factor Jo(2AQ5r/R), 
According to Eq. (8.73), the coolant temperature increases along the channel and 
reaches a maximum value 


71 WC 


(8.74) 


p 


at the exit of the channel This behavior of T h is illustrated in Fig. 8.8. The tem- 
perature of the fluid leaving channels other than the central hottest channel is also 
highest at the exit; but necessarily lower than the temperature of the coolant leaving 
the hottest channel. 

The temperature T c of the surface of the cladding can now be found as a 
function of position along the channel by observing that the heat transferred from 
a length dz of the fuel rod to the coolant is h€ c dz (T c - ?V), where C c is the cir- 
cumference of the clad rod. This is equal to the heat generated in the length dz of 
the fuel— namely, g^A jdz cos(xz/H). It follows that 

hCATc - T h ) = ql^Aj cos (^j . 
Inserting T h from Eq. (8.73) and solving for T c gives 

t c _ r M + h*?L [, + si „ f 5)1 + %lz cos («) . 

TtWCp L V///J hC € \ff J 



Figure 8.8 Axial temperature of 
fuel, T nu temperature of cladding 
surface, T c , and bulk temperature of 
coolant, as functions of distance 
along a coolant channel 
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This may also he written as 


TTWCp 


1 4- sin 


(!)• 


iiiiiir : 


where is the resistance for convective heat transfer (see Eq. [8.64]). 

The temperature T m along the center of the fuel can also he calculated 
Eq. (8.73). It is assumed that all of the heal produced in a slab of fuel of thick 
dz flows directly into the coolant. This rate of heat flow is equal to the temperature 
difference T m - T b divided by the total thermal resistance of a section of fuel and 
coolant dz thick. The resistance of such a section is simply the resistance calculated 
earlier for the full length of the fuel multiplied by the ratio dz/H. Equating h 
production and heat flux then gives 


H 


h) 


where R is given by Eq. (8,65) for plate-type fuel elements and by Eq, (8.67) or | 
Eq. (8.68) for cylindrical fuel rods. Solving for T m and introducing 7}> from Eq; 
(8-73) yields finally ~ |1F 


T m = 7k + 


71 WC 


P L 


1 + sin 


m 


+ <ltnm V f Rc ™ 


(f) 



A plot of Eqs. (8.75) and (8.76) is given in Fig, 8.8. It is evident from the 
figure that both T c and T m rise along the channel and reach maximum valuesin the 
upper pan of the channel beyond the midpoint of the fuel There are two reasons 
that the maximum temperature of the fuel occurs there, rather than at midchannci, 
where q m is the greatest. First, the temperature of the coolant continues to incre^jM 
past the midpoint. Second, the heat flux q" is determined only by the value of 
and this, being a cosine function, decreases very slowly in the vicinity of z ^ G§ 
But q", in turn, specifies the temperature difference T m - T b . Therefore, with 7}> 
increasing, T m must also increase to provide the actual value of q f \ Further along j 
the channel, q" begins to drop more rapidly, and T m eventually decreases. It is Uiis 
combined effect of a rising T h and a decreasing q" that gives rise to and determines 
the position of the maximum fuel (and cladding) temperature. 

To find the locations of these maximum temperatures, the derivatives of T c ou 
T rri are placed equal to zero. When this is done with Eq. (8.75), the result is easil 
found to be - • • - j f 

H 

Zc.m&s = ~CQC l (7tWC p R h ). 

n 



|:^fi|Jll 
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/Similarly, from Eq. (8,76), 


H 

Zmjmx = — COC* (XWC P R). 
71 
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(8,78) 


If The actual maximum values of T c and T m can be found by introducing Eqs, 
l| 77) and (8.78) into Eqs. (8.75) and (8,76). With a little trigonometry, it is easy 
" obtain the following results: 


I + 7 1 + a 1 


a 


where 


ahd 


*here 


a = nwCpRk 


7k max = Tio + flfmaxV/* 


i + yi 


0 


$ = TTWCpR. 


(8.79) 


(8.80) 


(8.81) 


(8.82) 


It must be emphasized that the previous calculations of 7i, F c , and T m as 
functions of 2 are highly approximate and give only the most qualitative behavior 
these temperatures. Furthermore, the results are not applicable to liquid metal 
■clolants. where there may be appreciable heat conduction in the z direction. The 
lerivations also do not apply to coolants that are undergoing bulk boiling, although 
t|ey are valid for locally boiling coolants. 

Example 8.7 

Pressurized water enters the core of the reactor described in Examples 83 and 8.5 
at a temperature of 543 D F and passes along the fuel rods at the rate of 3,148 Ih/hr 
per channel, (a) What is the exit temperature of the coolant from the hottest channel? 
(b) What is the maximum temperature of the cladding and the fuel in this channel? 

Solution. 


1. The volume of the fueled portion of the fuel rod is 

V f = 1.15 x I0" 2 ft 2 . 

From Example 8.3, = 4.66 x 10 ? Btu/hr-ft 3 , and from Table IV3, c p has 
a value of approximately 13 Btu/lb-°E Then using Eq. (8J4), 
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7ft.mn = 543 + — — t 1 

7T x 3148 x 1.13 
= 543 + 83 « 626°F. [i4«.l 

2. From Eq. (8.64) and Example 8.6, the total conveciive resistance is 

R h « J/7,500 x 2,T(0.01?5 + 0.002) x 12 = 9.07 x 10~ 5 T-hV/B|; , 

The location of the maximum cladding temperature is then given by Eq. (8.77): . 


= — cor'fr x 3148 x 1.3 x 9.07 x I0™ 5 ) 


SMI 


12 


= — cot _, (I.I7) = 2.71 ft. 


Thus, the maximum value of T c occurs 2.71 ft beyond the midpoint of i_ 
channel. The value of a in Eq. (8,80) is 1.17 and from Eq. (8.79), the maximum 
cladding temperature is 


Zv.nuK = 543 + 4.66 x 10 7 x 1.15 x 10~ a 

I + yT+(I.17F 
1. 17 

= 649°F [Mil* J 


x 9.07 x ur s x 


For the fuel, from Example 8.5, R f + R c = 6,17 x I0~\ so the total resistant I 
R is 6. 17 x 1 0^+9.07 x 10 5 = 6.26 x |0~ 3 °F-hr/Btu. Then using Eq&. (8.78) " 
(8.81), and (8.82). fi = 80.5, <,, Mtm = 0.047 ft (so that the temperature of t 
fuel is greatest at just about the midpoint of the fuel), and T m nm ^ 3,< 
\Ans.\ 

The Heat Transfer Coefficient— Nonrnetaliic Coolants 

The extent to which heat is transferred to a moving fluid, and hence the value 
of the heat transfer coefficient, depends on the details of the internal motions of 
the fluid as it flows along a coolant channel If every portion of the fluid moves 
parallel to the walls of the channel, then the heat travels radially into the fluid 
largely by conduction. The flow; in this case, is said to be laminar In contrast, 
if there are significant radial components of velocity fluctuations within the fluid, 
the heat is picked up at the wall by portions of the fluid and carried directly into 
the interior of the channel. This is the description of turbulent flow. Clearly, other 
things being equal, -heat is more readily transferred to a fluid undergoing turbulent ; 
flow than laminar (low. As a general rule, in those reactors in which the coolant i|| 
pumped through the system (as opposed to reactors cooled by natural convecUpn)| 
the coolant flows under turbulent conditions. 
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Figure 8.9 Temperature 
distributions in a nonmetallic 
coolant and in a liquid metal 
coolant undergoing turbulent flow. 


One of the consequences of the internal motions of the coolant undergoing 
turbulent flow is that the temperature distribution tends to be more or less uniform 
over much of the interior region of a coolant channel. Thus, as shown in Fig. 8.9, 
the temperature drops rapidly with distance in the vicinity of the fuel and quickly 
reaches the bulk temperature of the fluid. This is in marked contrast to the situation 
in conduction, where the temperature changes more uniformly with position. Even 
under turbulent conditions, however, some heat is always transferred to the interior 
of a coolant by conduction, but for the nonmetallic coolants this contribution is 
negligibly small. With ordinary water, for example, less than 1% of the heat transfer 
occurs by conduction when the water is undergoing turbulent flow. 

It is possible to characterize the flow of a fluid in terms of a dimensionless 
parameter known as the Reynolds number, which is defined as 


Re = 


D ( >vp 


(8.83) 


where D, is the equivalent diameter of the coolant channel discussed later, v is the 
average velocity of the fluid, p is its density, and f t is the fluid viscosity. The value 
of D,, is to be computed from the formula 

cross-sectional area of coolant channel 

(8.84) 


a = 4x 


wetted perimeter of coolant channel 

The term wetted perimeter in this expression refers to that portion of the perimeter 
of the channel that is structural and that therefore creates drag to passing coolant. 
For a hollow pipe carrying a coolant, the wetted perimeter is simply the interior 
perimeter of a section of the pipe perpendicular to its axis. Thus, with a circular 
pipe of inside radius a, the numerator in Eq. (8.84) is xa 2 and the denominator is 
2na, so that D e = 2a, the actual pipe diameter. For a bundle of rods of radius a in 
a square array of pitch s, as shown in Fig. 8.7, the cross-sectional area of a single 
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coolant channel is equal to s 2 - mr\ the wetted perimeter is lira, and Eqjflgjt 

It has been found experimentally that the flow of most fluids is laminar up 
to a value of the Reynolds number of about 2.000. Between 2,000 and 10,QOQ 4 " 
the How is partly laminar and partly turbulent, the fraction of each depending!; j| 
the structural details of the coolant channel the roughness of the channel wi jl 
and other factors. Above Re = 10.000. a fluid moves in fully developed turblHB 
flow. In summary, a high value of the Reynolds number implies a large amount of I 
turbulence, a high value of the heat transfer coefficient, and a high rate of heat flog 
into the coolant for a given difference in temperature between the cladding and the 
coolant. 

Example 8.8 | 

The fuel rods described in Example 8.5 are placed in a square array with a pitch 
of 0.600 in. The rods are cooled by pressurized water (P = 2,000 psi). which J 
flowing at a speed of 15.6 ft/sec. Calculate the Reynolds number for this coo^ 
flow assuming the water temperature is 60Q°E 

Solution. The radius of the fuel rods is 0.210 + 0.024 = 0.234 in. From Eq jf|t 
the equivalent diameter is 

, (0.6) 2 - ;r(0.234) 2 

O e = 2x 1 — ! _ —L _ 0,512 in = 0.0427 ft, 

n x 0,234 

The flow velocity is 1.5.6 x 3.600 = 56,200 ft/hr. From Table IV3, at 600 f F 
and 2,000 psi, p = 42.9 lb/ft 3 and fi = 0.212 Ib/hr-ft. Then substituting into Eq< 
(8.83) gives r§g§ 

0.0427 x 56. 200 x 42.9 
Re = _ « 486,000. 

The water is clearly flowing under turbulent conditions. 

;. iiiii: 

The numerical value of the heat transfer coefficient h, Eq, (8.60), is afuhcB0| w 
of the physical properties of the fluid, its rate of flow, and the diameter or effective; | 
diameter of the coolant channel These parameters, including /?, can conveniently 
be grouped together into three dimensionjess quantities — the Reynolds number, 
already discussed; the Nusselt number, Nu. defined as 
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(8.87) 


In these formulas, D e is the effective diameter of the channel, k and ju are the 
conductivity and viscosity of the fluid, and c p is the specific heat. All of these 
parameters must be specified, of course, in consistent units. 

It can be shown both from an analysis of experiments and from theory that 
convective heat transfer data can be correlated in terms of three dimensionless 
numbers — Re, Nu, and Pr. In particular, for ordinary water, heavy water, organic 
liquids, and most gases, all flowing through long straight channels under turbulent 
conditions, these data can be represented by an equation of the form 


Nu = CRe m Pr' ? , 


(8,88) 


where C, m, and n are constants. The value of h can be obtained from Eq. (8.88) 
using the definition of the Nusselt number Thus, 


Re^Pr*. 


(8.89) 


In using correlations of data of the type expressed by Eqs. (8.88) and (8.89)— 
or any correlations, for that matter — care must be exercised to determine the refer- 
ence temperature at which the fluid properties are to be evaluated. Frequently, this 
is the bulk temperature of the fluid 7*. It should also be noted that these equations 
are not valid for liquid metals, which must be considered separately. 

With ordinary water, heavy water, organic liquids, and gases flowing through 
long, straight, and circular tubes, the following values have been recommended 
for the constants appearing in Eqs, (8.88) and (8.89): C = 0.023, m = 0.8, and 
n = 0.4. [With the constants appearing in Eq. (8.88), Eq. (8.88) is known as the 
Dittus-Boelter equation,] These constants are often used in computing h for non- 
circular coolant channels by introducing the appropriate value of D e . However, 
significant errors in ft may result if the channels deviate substantially from the cir- 
cular configuration. 

In the important case of ordinary water flowing through a lattice of rods, 
parallel to the axis of the rods, recommended constants are: m = 0.8, n = 1/3, and 
C is given by 7 


C = 0.042- - 0.024 


(8.90) 


7 Webman, J., "Heat Transfer to Water Flowing Parallel to Tube Bundles," NucL ScL & Eng. 
6:78, 1979. 
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for square lattices with 1.1 < £ < 1 .3 and 

C = 0.026^ - 0.024 

for triangular lattices with U < £ < 1 .5. The quantities P and D are. respectively 
the lattice pitch and rod diameter. h 


Example 8.9 


-'illllll 


Calculate the heat transfer coefficient for the water flowing through the lauic^^B 
scribed in Examples 8.5 and 8.8. 

Solution, From Example 8.8, the lattice pitch is 0.6 in. the fuel rod radius is 0.234 

in. and C is Vf*" 

C = 0.042- - 0.024 = 0.0299. 

According to Table IV.3, c p = 1.45 Btu/lb-*F, // = 0.212 Ib/hr-fu anclf 
0,296 Btu/hr-ft^F. Introducing these values into Eq. (8.87) gives * ■ | ^ 

1.45 x 0.212 
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From Example 8.8. D e - 0.0427 ft and Re - 486,000. Equation (8.89) then g| Jjjj 

/i = 0.0299 ( ( ~~) (486. 000)°* 8 (l.039) ,/3 
= 7436 Btu/hr-ft 2 ~°R [/In.v ] 


The Heat Transfer Coefficient— Liquid Metals 

Heat transfer to liquid metal coolants is strikingly different from the transfef/1 
heat to ordinary fluids largely because the thermal conductivities of liquid meiaj| 
are so much higher than those of other types of coolants. At 400° E for instance, 
the thermal conductivity of liquid sodium is 46,4 Btu/hr-ft~°R whereas it is on|§ 
0.381 Btu/hr-ft-°F for ordinary water and 0.1 15 Btu/hr-ft- 0 F for helium at 1 atriti||§f 
Thus, at this temperature, the conductivity of sodium is 122 times that of water anf| 
400 times that of helium. 

One important effect of the high conductivity of liquid metals is that, 
when they are flowing under turbulent conditions, these coolants absorb Hcial 
mostly by conduction. This is in marked contrast to the situation with nonmeta|j 
lie coolants, where, it is recalled, heat transfer occurs largely as the result of *M 
internal motions of the fluid. 
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i Since heat flows into a liquid metal primarily by conduction, the temperature 

I distribution within a coolant channel containing a liquid metal resembles the tern- 
fcy perature distribution in a solid conductor whose axis and circumference are held 
at different temperatures. Thus, as shown in Fig. 8.9, the temperature varies more 
slowly across the channel with a liquid metal than with a nonmetallic coolant. In 
I both cases, the bulk temperature of the coolant is given by Eq. (8.6 1). 

The heat transfer coefficient for liquid metals has been studied extensively 
in recent years, and correlations have been given that cover most of the situations 
4 encountered in practice. For the case of a liquid metal flowing under turbulent 
|; conditions through a hexagonal lattice of rods, parallel to the rods, Dwyer (see 
§ References at the end of the chapter) has given the following correlation: 

§ Nu = 6.66 + 3 J 26($/d) + 1 AMts/d) 2 + 0.01 55(*Pte) a86 f (8.92) 

where s/d is the ratio of lattice pitch to rod diameter, * is a function given graph- 
( ically by Dwyer that has been fitted by Hubbard 8 by the expression 

I ~r , 0.942(5 /rf) J 4 

I PctRe/lO 3 )'- 2 * 1 ' { } 

and Pe is the so-called Peclet number, which is given by 

Pe = Re x Pr = D ' vpc P , (8.94) 

A. 

Equation (8.92) is valid only for lattices with s/d > 1 .35. Values of the Nusselt 
number for more closely spaced hexagonal lattices are given in tabular form by 
Dwyer, Berry, and Hlavac (see References). 

Equations (8.92) and (8.93) may be used for square lattices by replacing the 
ratio s/d in these expressions by 1 .075(f /rf)„ where (s/d) s is the value of s/d for 
the square array. For rightly packed square lattices, the following correlation may 
be used: 

Nu = 0.48 + 0.0 1 33(Pe) a7 °. (8.95) 


8.5 BOILING HEAT TRANSFER 

Up to this point, it has been assumed that a liquid coolant does not undergo a 
change in phase as it moves along a coolant channel absorbing heat from the fuel. 
It was assumed that the coolant does not boil. However, there are some distinct 
advantages in permitting a reactor coolant to boil. For one thing, the coolant pres- 
sure is much lower when the coolant is allowed to boil than when boiling must be 

a E R. Hubbard 111, private communication, 
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prevented. In addition, for a given flow rate and heat flux, lower cladding and fuel 
temperature are required for a boiling than for a nonhoiiing coolant. For these real 
sons, boiling of a restricted nature is now permitted in pressurized water :mactorsf iff 
although steam is not produced directly in these reactors. With boiling water reai< 
tors not only is advantage taken of higher heat transfer rates, but, as pointed out | 
Chapter 4, by producing steam within these reactors, the entire secondary coolant 
loop of the PWR can be eliminated 

Boiling coolants other than ordinary water have been considered in a 'nunf - 
her of reactor concepts, but none of these has reached a practical stage of deyelop|^r 
menu Therefore, the following discussion pertains largely to water-cooled reactors, 
However, the principles to be considered also apply to other types of liquid-cooied 
reactors. 

Boiling Regimes 

To understand the phenomenon of boiling heat transfer, consider an experiment/ 
which the heat flux (Btu/hr-ft 2 or W/cnr ) from heated fuel rods to a flowing liqui 
coolant is measured as a function of the temperature of the surface of the rods for a 
given system pressure and flow rate. The results of such an experiment are shown 
in Fig. 8.10. As indicated in the figure, the heat flux increases slowly as the rod 
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Figure 8.10 The logarithm of the heat flux into a flowing coolant as a 
function of the surface temperature of the coolant channel. 
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|r temperature is increased at low values. In this temperature range, between points A 
|j: and B in Fig. 8.10, heat is transferred to the coolant by ordinary convection with no 
|j change in phase, and the heat transfer coefficient is determined by the correlations 
I 'given in the preceding section. 

; |lv As the surface temperature of the fuel is increased further, a point is eventu- 
ally reached where bubbles of vapor begin to form at various imperfections on the 
| surface of the fuel rods. This occurs at about the point B in Fig. 8. 10, and is a form 
tS of boiling called nucleate boiling. As the bubbles are formed, they are carried away 
It-from the rods and into the body of the coolant as a result of the turbulent motions 
I lof the fluid. However, as long as the bulk temperature of the coolant is less than its 
i ; saturation temperature, the vapor in the bubbles soon condenses to the liquid state 
4 and the bubbles disappear from the coolant. There is no net production of steam 
y under these circumstances, and the boiling process is termed subcooled nucleate 
%. boiling or local boiling. If and when the bulk temperature of the coolant reaches its 
saturation temperature, the bubbles persist within the coolant stream, there begins 
* to be a net output of steam, and the system is now said to be undergoing saturated 
f nucleate boiling or hulk boiling, 

The onset of local boiling and the transition to bulk boiling are shown in the 
lower two portions of Fig. 8. 1 1 for a fluid flowing vertically through a heated pipe. 
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Figure 8*1 1 Flow patterns in a 
vertical heated channel (From L. S. 
Jong, Boiling Heat Transfer and 
Two-Phase Flow. New York: Wiley, 
1967.) 
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With local boiling, as indicated in the figure, bubbles exist only near the surface of 
the pipe, and most of the pipe is filled with liquid. In the region with bulk boiling, 
however, the bubbles are distributed throughout the fluid, which is said to move in 
bubbly flow. Under certain conditions— namely, at high-flow velocities and a large 
concentration of bubbles (large void fractions)— the bubbles combine to forma 
void space along the center of the channel. This is called annular flow. 

In any event, with the onset of nucleate boiling, heat moves readily into the 
liquid. At every temperature in this region between B and C in Fig. 810, heat 
transfer is more efficient than by ordinary convection. There are two i^on^ljfte 
this. First, heat is removed from the rods both as heat of vaporization and as sen- •; 
sible heat. Second, the motions of the bubbles lead to rapid mixing of the fluid. 
The rapid increase in the heat flux with temperature, shown by the steep slope of 
the nucleate boiling region in Fig. 8,10, is explained by the fact that the density 
of bubbles forming at and departing from the rod surface increases rapidly with 
surface temperature . . || 

For ordinary water at pressures between 500 and 2,000 psia undergoing "nucle 
ate boiling, either local or bulk, the heat flux can be computed from the follov 
correlation: 9 


60(?7 10 6 ) 1/4 



where T c is the surface temperature of the cladding (°F). T m is the saturation tern* 
perature (°F), P is the system pressure (psia), and q' is the heat flux in Btu/hr-fr. 
Equation (8,96) is known as the Jens and Lottes correlation and is valid 
channel geometry. This correlation is not valid, however, if the coolant is tm< 
ing annular flow. 

It is important to estimate the fuel surface temperature T LB at which Iocs 
boiling begins to know which correlation to use for calculating heat flow into " 
coolant— a correlation for convective heat transfer or the Jens and Lottes equati 
Although there is no single, fixed temperature at which boiling starts, 71b ts-ii' 
computed from the simple formula 

Tub sb T sm 4- (T € — T m ht - — . 

h 

In this equation, T %m is again the saturation temperature, (T c - T m )st is -; the 
ference between the cladding and saturation temperatures as computed ^!|f|§ 
Jens and Lottes correlation, q h is the heat flux, and h is the convective heat fifcnsfe* 
coefficient, 'f^ 

The application of some of the prior results is illustrated in Fig. 8.12, wJi&h 
shows the parameters of the water coolant in a boiling-water reactor as a functtp* 1 * 

9 W. H. Jens and P. A. tones, USAEC Report ANL-4627, 1951. 
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Figure 8.12 Thermal analysis of the hottest channel of a BWR. Here, 
x refers to void fraction not quality. (Courtesy of F. R. Hubbard ILL) 


of distance up the hottest channel In this reactor, the fuel rods are about 12 ft long 
and the system pressure is 1,035 psia. At this pressure, the saturation temperature 
of water is 548 .7°F, which can be estimated from the data in Table IV.2. The lower 
curve shows the acmal heat flux into the channel. The skewing of q" away from 
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the bottom of the core is due to the presence of control rods in this region The cal 1 
ciliated value of T XM and the outer cladding temperature / . are (he saine-namelv 
5^8 F approximately 13 in up the channel. Local boiline commences at this po j m ' 
The bulk temperature of the water increases stead.lv from its entering valu|§f 
526'F at the bottom of the core and reaches the saturation temperature of MS ?F 
about 45 in along the channel. It is at this point that bulk boiline begins Thus th- 
heat transfer ,s convectivc for about the first 13 in up the channel, local boiling 
occurs between 13 in and 45 in. and throughout the rest of the channel the water 
undergoes bulk boiling. (The flow never becomes annular.) 

It should be observed in Fig. 8. 1 2 that the difference in temperature between 
the cladding and the coolant, after rising in the convective region in response to 1 
increasing q , drops sharply after the onset of nucleate boiline at 13 in, allhou* 
q is still increasing. This illustrates how readily heat flows into a boiling coolam 

It should also be noted in Fig. 8.12 that the void fraction of the coolantlflj 
to rise from zero at about the point where local boiline begins, although iherfii 
no net production of steam until the onset of bulk boiling further up the channel 
The reason tor this is that, although the steam bubbles thai form at the rod surface 
collapse as they pass into the interior of the channel, their place at the surface is 
immediately taken by other bubbles. As a consequence, there is alwavs a steady 
state distribution of voids across the channel. 

The Boiling Crisis 


•a 

■■I 


llii! 


To return to the discussion of the experiment depicted in Fig. 8.10. if the temper- 
ature of the fuel surface is increased in the nucleate boiling region, the density of ' 
bubbles at or near the surface of the fuel rods also increases. Eventually however, 
a point is reached where the bubble density becomes so great that adjacent bub- 
bles coalesce and begin to form a vapor film across the surface of the rods At this 
point, which corresponds to the point C in Fig. 8.10. the .svstem is said to he in a 
boilmg casts or in a condition leading to a departure fmm nucleate hailing (DNB). 
The heat flux at or just before the boiling crisis is called the critical heal flux- 1 
(CHF) and denoted as q c . Sometimes, for reasons that are evident momentarily, th- 
is called the burnout flwe. 

With the onset of the boiling crisis, the heat flux into the coolant begins to 
drop. This is due to the fact, that over the regions of the rods covered by sapor 
film, the heat is forced to pass through the vapor into the coolant bv conduction 
and radiation, both of which are comparatively inefficient mechanisms for heal 
transfer. The heat flux continues to drop, more or less erratically (as indicated by? 
the dotted lines in Fig. 8.10), with increasing fuel temperatures as the total arelbl 
the him covering the feel increases. In this region of Fig. 8.10. the svstem is said 
to be experiencing partial film boiling. 
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Eventually, when the rod surface temperature is high enough, the vapor film 
|B| )Vcrs the entire rod and the heat flux to the coolant falls to a minimum value (point 
Beyond this point, any increase in temperature leads to an increase in the heat 
Ipflflux simply because heat transfer through the film, although a poor and inefficient 
% process, nevertheless increases with the temperature difference across the film. The 
svstem, in this case, is said to be undergoing full film boiling. 
* The existence of the various boiling regimes and the boiling crisis is an im- 
portant consideration in the design of a liquid-cooled reactor, as may be seen from 
the following example. Suppose that a water-cooled reactor is designed and/or op- 
erated in such a manner that at some point along a coolant channel the water un- 
dergoes nucleate boiling near DNB conditions. If the reactor power is suddenly 
increased so that the heat flux into the water rises above the DNB value of q c> 
partial film boiling will immediately begin in this channel However, as explained 
previously, the formation of the film impedes the transfer of heat to the coolant. 
i :As a consequence, the heat confined, so to speak, within the fuel raises the fuel 
temperature and the surface temperature of the rods — forming the channel. This, 
iit^rn, leads to an increase in the area of the film, which leads to a further de- 
ft- crease in the heat flux, a further increase in rod surface temperature, and so on. In 
this Way, the wall temperature rapidly increases along the boiling curve from point 
C toSpoint E, Long before £ is reached, however, the temperature of the fuel will 
taiam-such high values (several thousands of degrees Fahrenheit) that the fuel will 
c /partially melt, the cladding will rupture, and fission products will be released into 
fflpioJant... As noted earlier, these are occurrences that must be prevented at all 
costl? For this reason, it is important to know the value of q c and to keep a reactor 
m operating near the DNB point. 

It must not be implied from Fig. 8.10 that it is necessary to have saturated 
nucleate boiling (bulk boiling) before the onset of the boiling crisis. The figure was 
Tfeten'ded merely to show that there are two types of nucleate boiling and that it is 
possible to make the transition from subcooled to saturated boiling if the bulk tem- 
||j|||f^ire;of the coolant exceeds the saturation temperature. Bubbles at the surface 
of r the fuel can combine to form areas of vapor film, the beginning of the boiling 
crisis, although the bubbles would condense if they traveled to the interior of the 
.coolant stream. 

/ |Jq|\ great many correlations have been developed from data on q c {see long in 
ijMep^ierences for a tabulation of these correlations). As would be expected, these 
||eor|elations differ depending on whether the boiling crisis is reached from sub- 
* §cooled or bulk boding conditions. For subcooled boiling, the following correlation 
,.yby Jens and Lottes 10 has been widely used: 


fW. H. Jens ami P. A, Lottes. ibid 


_ 




If 
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TABLE 8.1 PARAMETERS FOR JENS AND 
LOTTES CORRELATION 


P (psia) 

C 

m 

500 

0.817 

0.160 

1,000 

0.626 

0.275 

2,000 

0.445 

0.500 


q'' = Cx 10 6 


(lO 6 ) 


at; 


•0.22 
sub > 


where C and m are pressure-dependent parameters given in Table 8.'lM 
coolant mass flux in lb/hr-ft 2 , and Ar sub is the difference between toetsafssa 
and local temperatures in °F. 

Another correlation for subcooled boiling was obtained by BematKSS? 
a combination of the three equations: 


q"c = h c (T wc - T b ) 
T wc = 102.6 In P- 


97.2P 
P + 15 


- 0.45v + 32 


hc = 10,890 • 


IP 

% 


(81 


In these expressions, T wc is the wall (cladding) temperature at the onset^|L_ 
boiling crisis, T b is the bulk temperature, P is the pressure in psia, v is the coolant 
velocity in ft/sec, D e is the equivalent diameter in feet, and A is defined^ 
heated perimeter of a channel in feet divided by n. The Bernath correlation is||j 
for pressures between 23 and 3,000 psia, fluid velocities between 4.0 and 54./ 
and for D e between 0.143 and 0.66 in. 

Example 8.10 

For the PWR lattice described in Examples 8.3 through 8.8, calculate for the onset 
of the boiling crisis the: (a) cladding temperature, (b) heat transfer coefficient, "and 
.(c) critical heat flux. iff 


Solution. 

1. Using Eq. (8.99) of the Bernath correlation, 

97.2 x 2,000 


if 


T wc = 102.6 In 2,000- 


2,000+ 15 
n Bernath, L., Transactions A ./. Ch.&, 1955. 


- 0.45 x 15.6 + 32 = 708°F. [AnsHtfk 
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2. From Example 8.8, D e = 0.0427 ft. The heated perimeter of the channel is 
2?ra, where a = 0.234 in. Thus, 


D, ; = 2x0.234/ 12 = 0.039 ft. 

Then from Eq. (8.101), 

h c = 10,890 x 


0.0427 


0.0427 + 0.039 + (0.0427) 06 
= 10,659 Btu/hr-ft 2 -°F. [Ans.] 
3. With a bulk water temperature of about 600°F from Example 8.6, the critical 
heat flux is 

q l = 10, 659 x (708 - 600) = 1.15 x 10 6 Btu/hr-ft^F. [An,] 

E q uations(8.89)tr^^ 

resents the mass of ^^^^^^^^Ja the case of the 
more or less through the center of & General Electric 

approach to DNB from bulk boiling conditions workers at tn ^ 

Company (GE) ^*&t£^£££Z de/end on the^ 

mass flow rate of vapor (lb/hr-ft 2 ) (8 . 10 2) 
* = "mass flow rate of vapor-liquid mixture (lb/hr-ft 2 ) 
The GE correlation, due to Janssen and Levy, 12 is as follows: 

iL = 0 705 + 0.237(G/10 6 ), X < Xi 
10 6 

£- = 1 634 - 0.270(G/10 6 ) - 4.710 X , Xi < X < X2 
10 6 

iL = 0 605 - 0.164(G/10 6 )-0.653 X , X2 < X, 
10 6 


48 x 15.6 


where 


Xl = 0.197 - 0.108(G/10 6 ), 
X2 = 0.254 - 0.026(G/10 6 ), 


(8.103) 
(8.104) 
(8.105) 

(8.106) 
(8.107) 


Accn W&l 1962. The General 
HE. Janssen and S. Levy, General J^^SSSto. due to Hench and 

Electric Company has also published a poruon of a more recent (1966) 

Levy. 
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and G is again the mass flux in lb/hr-ft 2 . These equations refer to a system pressure! 
of 1 ,000 psia. For other pressures, q" c can be found from ' || 

q" c (P) = q" c (1,000 psia) + 400(1,000 - P). (8.108§i 

The GE correlation is valid for pressures between 600 and 1450 psia, G from 0 4 M 
10 6 to 6.0 x 10 6 lb/hr-ft 2 , values of X up to 0.45, D e from 0.245 to 1.25 in, and foil 
coolant channels between 29 and 108 in. in length. U 
All of the foregoing correlations for the critical heat flux were developed! S 
largely from data for circular pipes. When used to determine q" c for noncirculafe 
channels or for coolant flowing through lattices of fuel rods, these correlations 
can be counted on to provide only the roughest approximations for q" c . Various 
correction factors for the various geometries have been developed, however, buf • 
these are too detailed to be discussed here (see Tong in the References). 

8.6 THERMAL DESIGN OF A REACTOR 

As pointed out earlier in this chapter, reactors must be designed in such a way thafe^ 
the fission products remain confined within the fuel at all times— throughout the f 
operational lifetime of the core, during shutdown, and under accident conditions!^ 
when the fuel may be denied normal cooling. This places upper limits on the VemM 
perature of the fuel and/or its cladding. The currently accepted design criterion is% 
that the integrity of the cladding must be maintained through all operating con|| 
ditions. Since the expansion of the fuel on melting can rupture the cladding, tWs\ 
design criterion is essentially equivalent to the requirement that the fuel must notl ' 
melt. P 

The melting point of U0 2 which depends somewhat on burnup— that is, on^ 
the fraction of the fissile atoms that have undergone fission— is generally between^ 
5,000°F and 5,100°F (2,760 °C and 2,815°C) for commercially available U0 2 . hf" 
most reactors fueled with U0 2 , the maximum fuel temperature is somewhat below!; :k 
4,500°F (2,480°C). " 

In the HTGR, where the fuel consists of small particles of uranium and tho|$ 
rium dicarbide in a carbonaceous binder, the maximum permitted fuel temperature^., 
is about 6,500°F (3,600°C). * p 

Natural or enriched uranium metal melts at 2,070°F (1,132°C), but it also^> 
undergoes two changes in solid phase (alteration in crystalline structure)— first t| 
at 1,234°F (668°C) and again at 1,425°F (774°C), However, above about 750°%? 
(400°C), the strength of the metal decreases rapidly. This permits fission product^ 
gases to gather and diffuse into pockets in the fuel, which in turn leads to the ex|* 
pansion of the fuel and to cladding failure. Therefore, it is usual practice to design^ 
uranium metal-fueled reactors with maximum fuel temperatures below 750°F. This§| 
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low temperature is not as serious a drawback in transferring heat to the coolant as 
it would at first appear since the thermal conductivity of uranium metal is much 
higher than, say, that of U0 2 . At 600°F, for example, k f is 6.5 times higher for 
uranium than for UO2. 

The DNB Ratio 



To prevent the penetration of the cladding at any point in a water-cooled reactor 
due to the onset of film boiling, this type of reactor must be designed so that the 
heat flux q" is always below the critical (burnout) heat flux q c . For this purpose, it 
is convenient to define the DNB ratio as 


DNBR = 


factual 


(8.109) 


In this ratio, q c is the critical heat flux computed as a function of distance along 
the hottest coolant channel from the appropriate correlation given in the preceding 
section, and # actual is the actual surface heat flux at the same position along this 
channel. If heat flow parallel to the fuel can be ignored, then 


*2 actual q > 


(8.110) 


where q f " is the heat production per unit volume, Ay is the cross-sectional area of 
the fuel, and C c is the circumference of a heated rod. These two functions and the 
DNB ratio are shown in Fig. 8. 1 3 for a 1 ,500-MW boiling- water reactor. BWRs are 
currently designed with a minimum DNBR of 1.9. 13 By setting this lower limit on 
the DNB ratio, there can be reasonable certainty that burnout conditions will never 
be reached anywhere in the reactor, even during transient, overpower situations. 
The minimum design value for the DNBR for PWRs is 1.3. Needless to say, the 
establishment of a minimum DNB ratio provides a major limitation on the design 
of water-cooled reactors. 

Hot Channel, Hot Spot Factors 

With water-cooled reactors, the maximum heat flux anywhere in the core is limited 
by the DNB ratio. In gas-cooled reactors, the maximum q " is determined solely by 
the requirement that the fuel temperature remain well below the melting point for 
this value of q^. In any reactor, the extent to which q' max exceeds the average heat 
flux in the core is given by the hot channel factor, also called the hot spot factor 


13 The minimum DNB ratio is also called the minimum critical heat flux ratio (MCHFR). 
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Figure 8.13 The actual heat flux, the computed critical heat flux, and 
the DNB ratio for a boiling-water reactor. (Courtesy of General Electric 
Company.) 


This is defined by the relation 


F = 


where q^ is the average heat flux in the core. 

There are several reasons that differs from and therefore why |gp 
different from unity. The most important of these is that, as seen earlier, the pojjKe§ 
distribution across the core is not flat. If this were the only reason, F coul|^ 
computed directly for any given reactor design in more or less the same wa||§£ 
the maximum-to-average flux ratio Q was calculated in Chapter 6. Thus, ?LJ|p 
function of various nuclear parameters that determine the power distribution,^^ 
q[ v is given by the formula 

ql = Pj A, 

where P is the thermal power output of the core and A is the total heat transfer^ 
of the fuel (with cladding). The previous computation yields what is known afflfi 
nuclear hot channel factor F N . U " ' 

14 For a uniformly fueled reactor at the beginning of core life and in the absence of cohjgl 
rods, F N is the same as Q. 
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In addition to the effect of the nonflat power distribution, q^ may differ from 
the computed value of q^ as the result of various statistical factors over which 
the reactor designer has litde or no control. For example, the amount of fissile 
material included in the fuel pellets of a BWR or PWR at the time of manufacture 
varies slightly from pellet to pellet because of the inherendy statistical nature of 
the manufacturing process. Pellets containing more fissile material produce more 
power. If such a pellet was located at the point where q" is highest, the value of 
9max wou ld be higher than calculated. 

In a similar way, manufacturing tolerances in fuel assemblies may result in 
slight bowing of the fuel rods, leading to reduced coolant flow and excessive heat- 
ing of a portion of the fuel. Similarly, fluctuations in the thickness of the cladding 
may give rise to hot spots where the cladding is thinnest. 15 Also, certain aspects of 
the coolant flow are inherently statistical in character and tend to give fluctuations 
in the heat flux. 

These various mechanisms, and others by which q mz ^ can differ from its com- 
puted value, taken together are described by the engineering hot channel factor F E . 
The overall hot channel is then 


Each of the individual mechanisms is described, in turn, by an engineering hot 
channel factor, Fe, x - 

The engineering subfactors are obtained from data on fabricated reactor com- 
ponents or from tests, such as coolant flow tests on assembled portions or mockups 
of the reactor. Consider, for instance, the amount of fissile material mi per unit 
length of the fuel rods. For a given thermal flux, q^ is very nearly proportional 
to m/. When measurements of m/ are carried out on manufactured rods, it is found 
that the actual values of mi form a normal distribution about some mean or average 
value nil like that shown in Fig. 8.14. With such a distribution, all values of mi are 
presumably possible. However, specific values of m/ become less and less likely 
with increasing deviation from the mean. In particular, it can be shown that the 
probability that m/ exceeds m x by more than 3a, where a is the standard deviation 
of the distribution, is only 1.35 per 1,000 — an unlikely occurrence. 

In those situations like the one just described, in which q^ is proportional to 
a normally distributed engineering variable x, the engineering subf actor is defined 
as 


F — F^Fe- 


(8.113) 


Fe, x = 


x + 3cr(x) 



(8.114) 


x 


15 This is because thin cladding leads to a larger gap (and hence a larger AT between the fuel 
pellet and the cladding). 
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Figure 8.14 Normal distribution of measured values of w/. 


Here x is the mean value of x and a(x) is the standard deviation in the measiim| 
ment of x. For example, in the case of the amount of contained fissile material; tHe 
subfactor is fe$ 


FE,mi = 1 + 


3a (m/) 


Example 8.11 

In a certain PWR, the average linear density mi of U0 2 in a fuel rod is 0.457 lb/ft? 
The standard deviation in the measured values of m/ is 0.0122 lb/ft. What is F E Jtl 

Solution. It is merely necessary to substitute these values into Eq. (8. 1 15). This gives 


Fe,t7h = 1 + 


3 x 0.0122 


0.457 
= 1.08. [Arts.] 


This result means that, due to statistical fluctuations in the amount of fissile material! 

included in fuel rods during manufacture, there are 1.35 chances in 1,000 that 

will exceed its computed value by more than 8%. Jp 

Calculations of other engineering subfactors, especially those that do not in|| 
volve a direct proportionality between q' max and the statistical variable, are more|f 
complicated and require an analysis of the relationship between the variable in; 
question and # max . As noted earlier, such subfactors are often obtained from engif 
neering test data. 

Once the various engineering subfactors have been determined, it is possible^ 
to calculate F £ . This can be done in a number of ways. The most obvious thing is to| 


fife* 
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•f t u itiply all of the subfactors together, and this was done in the case of some early 
W i actors. However, the resulting value of F E leads to an unnecessarily conservative 


W acto f design since, in using the product of subfactors, it is implicitly assumed that 
|jf eV eral statistically unlikely events occur simultaneously. The preferred method for 


(jiputing is based on an analysis of overlapping statistical distributions that 
f|^ to0 complex to be reproduced here (see Tong and Weisman in the References). 


<^ e edless to say, with F Ey 


known, F E can be computed. 


i 


| pe actor Design 


designing a reactor, there is considerable overlap between those aspects of the 


f0 blem that are specifically nuclear — the neutron flux and power distributions, 
ent behavior, and so on — and those factors involving the removal of heat 
the system. The final reactor necessarily represents a series of compromises 


%etw een nuc l ear thermal-hydraulic factors, all made in the interest of safety 
|Lfld economy. One noted reactor designer frequently refers to reactor design as 
fpthe art of compromise." 16 

% Of course, most of the reactors manufactured today are not designed from 
it cfa tch. They are either identical to or improved versions of existing plants. Only 
finely is an engineer called on to design a reactor from the beginning. Nevertheless, 
!f the several steps in such a design study are now considered to illustrate some of the 
| c0 ncepts developed earlier. 

% Consider the design of a specified type of power reactor having a given ther- 
f 0ial poW er output. The first problem is the selection of the core materials. These 
f ^ e usually chosen on the basis of previous experience, availability, and cost. The 
ll amount and enrichment of the fuel are determined in a preliminary way by various 
$ fl uclear considerations, including estimates of core lifetime and radiation damage 
if -jo the fuel elements. 

With this preliminary choice of materials, a fuel lattice arrangement is next 
a proposed. The operating parameters of the reactor — the neutron flux, the power and 
4 t6 mperature distributions, the coolant density and/or void distribution, the temper- 
Jl ature and heat flux along the hot channel — are then computed for the proposed 
?i lattice. j t should be noted that such computations are inherently iterative, at least 
If for some types of reactors. With the PWR or BWR, for instance, it is necessary 
€ to know the water density as a function of position to compute the neutron flux 
3$ and, from this, the power distribution. However, the power distribution must be 
i known before the water density can be found since the density is determined by 
the heat flux distribution, which in turn depends on the power distribution. What 


l6 B. Minkler, private communication. 
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must be done in this case is to assume some density distribution to start with4%\ 
unity throughout the core. The flux and power distributions are then computed^ 
this water density, and a new density distribution is determined from the com$l§ 
heat flux. The neutron flux and power distributions are now computed usinjfif 
new density distribution, and this in turn leads to a third density distribution^ Tfle 
calculations are continued in this way until convergence is obtained— that is,^ffii 
the density distribution provides a power distribution that gives the same denl| 
distribution. Computer programs are available by which such iterative calculatifs 
can be carried out automatically. Incidentally, during this preliminary stage in the 
design of the core, the neutron flux is usually obtained from a multigroup diffu- 
sion calculation using only a few groups. To reduce computation costs, multig^ 
calculations with many groups are normally performed in the final stages ofM 
design. ~ ff 

If the reactor is water cooled, the DNB ratio is next computed as a functidnfl 
distance along the hottest channel. As pointed out earlier in this section, the mini- 
mum DNBR must be equal to some preset limiting value. If the computed DNBR 
is not equal to this value— it would only be luck if it were at this point— then|t§ 
parameters of the lattice or the coolant flow are changed, the previous calculation! 
of the core are repeated, and a new DNBR is obtained. These computations S| 
repeated until the desired (DNBR)^ is found. This also determines the maximum 
heat flux, although, as indicated in Fig. 8.13, q' max does not necessarily occur k 
the point where the DNBR is smallest. If the reactor is not water copied, q^j^s 
determined solely by the limiting value of the fuel temperature. For either type$f 
reactor, water cooled or not, a preliminary calculated value of # max has now been 
found. ^ " >l " K |f 

At this point, the overall hot channel factor is obtained using the methods 
discussed earlier. The average heat flux throughout the core can then be computet 
from Eq. (8.111): 


4av = 


"max 


i 

(8.1 li- 


ft.? 

By the simple device of designing the reactor to operate at this average heat fluff 
the actual value of q' max is reduced below its calculated value by the factor F E . Thi| 
provides the reactor with the desired margin of safety. || 
Once q" w is known, the total heat transfer area can be found from Eq. (8. 1 12'jjf 
Finally, from the heat transfer and the dimensions of the fuel rods, it is possible to? 
compute the total number of rods required in the reactor. This number will undoubfe 
edly be different from the one used at the start of the design, and therefore many of£ 
the foregoing calculations must be repeated until a completely self-consistent re|P 
actor is obtained. At this point, detailed multigroup and two-dimensional neutron; 
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diffusion calculations can be undertaken, along with various transient analyses and 
lifetime-burnup computations, to refine the design. 

Example 8.12 

The core of a PWR consists of a lattice of fuel rods 12 ft long and 0.5 in. in diameter. 
The reactor operates at a thermal power of 3,000 MW; it has a maximum calculated 
heat flux of 539,000 Btu/hr-ft 2 and an overall hot channel factor of 2.80. Calculate 
(a) the total heat transfer area; (b) the number of fuel rods in the core. 

Solution. 

1. The average heat flux is given by Eq. (8. 1 12): 

539,000 


2.8 


= 193,000 Btu/hr-ft 2 . 


With a reactor power of 3,000 MW = 3,000 x 3.412 x 10 6 Btu/hr, the heat 
transfer area is 

3,000 x 3.412 x 10 6 9 
A = mM = 53,000 ftMAn,] 

2. The surface area of each fuel rod is 12 x (0.57r/12) = 0.57T ft 2 . The total 
number of rods is then 

53,000 

n = — = 33,740. [Ans.] 
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